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Abstract 



^ We present a new classification of elementary cellular automata. It 

2 is based on the structure of the network of states, connected with the 

^ transitions between them; the latter are determined by the automa- 

—L ton rule. Recently an algorithm has been proposed to compress the 

^ network of states (M. J. Krawczyk, Physica A 390 (2011) 2181). In 

Q this algorithm, states are grouped into classes, according to the local 

O symmetry of the network. In the new classification, an automaton is 

described by the number of classes #(iV) as dependent on the sys- 

1—1 tern size TV. In most cases, the results reflect the known classification 

^ into 88 groups. However, the function #(iV) also appears to be the 

same for some rules which have not been grouped together yet. In this 

way, the automaton 23 is equivalent to 232, 77 to 178, 105 to 150, the 

l^ pair (43, 113) to the pair (142, 212) and the group (12, 68, 207, 221) to 

• the group (34,48,187,243). Furthermore, automata 51, 204, the pair 

^ij (15, 85) and the pair (170, 240) are all mutually equivalent. Results are 

(^ also presented on the structure of networks of states. 



1 Introduction 

Attempts of cellular automata classification were made by several authors. 
Proposed approaches are based, either on observations of the automata be- 
haviour, or properties of the rules themselves. The most famous classifica- 
tion of the elementary cellular automata CA belongs to the first category, 
and was introduced by Wolfram [1^ . The classification regarded the character 



of automata evolution and led to the division of the whole set of elemen- 
tary rules into four classes. The first group includes rules which applied to 
a given sequence, cause removal of any randomness, and as a consequence 
lead to a stable state. The rules included into the second class lead to a sta- 
ble or oscillating character of evolution. The third class covers rules which 
show with pseudo-random or chaotic behaviour. The fourth class is espe- 
cially interesting as it covers rules which lead to complex evolution. The 
classification makes it possible to see differences in behaviour of the system 
observed during the evolution ruled by the particular automata. 

Besides the Wolfram's classification there are also different propositions. 
One of them was introduced by Li at al. [2] . The classification is based on 
the observation of the asymptotic behaviour of the elementary CA, and, as 
a result, six types are distinguished. CA for which evolution leads to: a) 
spatially homogeneous fixed point, b) spatially inhomogeneous fixed point 
or a uniform shift of a fixed pattern, c) periodic behaviour or shifted periodic 
behaviour, d) locally chaotic behaviour, e) chaotic behaviour or f) complex 
behaviour. Another classification is based on the concept of equicontinuity 
which originally was introduced by Oilman [3] and then extended by Kurka 
t4j , where an equicontinuity set contains a set of its Lyapunov stable points. 
In this scheme, similarly to Wolfram's classification, each OA can be in- 
cluded into one of the four categories: equicontinuous OA, OA with some 
equicontinuous points, sensitive but not positively expansive and positively 
expansive. It is not clear a priori to which category a particular automaton 
belongs, but for the elementary OA this problem can be solved [5]. A re- 
view of the classification methods, along with a discussion of problems that 
arise when attempting classification of automata, is presented by Sutner in 
[6]. Analysis of cellular automata properties is an important problem as 
they have a lot of applications in modelling of physical processes, such as 
diffusion, phase transition, wave propagations or road traffic analysis {7}. 

In this paper we propose a method of class identification which is based 
on the similarities observed in patterns of connections between possible 
states of the analysed system. Here we say that two states are connected 
if one state is obtained from another one, according to a given automaton 
rule. The presented method takes into account the symmetry observed in 
the state space, and it allows us to indicate equivalent - according to symme- 
try - rules. The method is exhaustive, and is based on scrupulous analysis 
performed for sequences of lengths up to A^ = 19 (in some cases longer if 
the character of results does not allow conclusions to be drawn earlier). 
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Table 1: Rule equivalence. The bar means negation, i.e. = 1 and 1 = 0. 

The paper is organised as follows. In the next section the analysed 
system is presented. In the section after that our procedure for identification 
of classes is described. Section |4]is devoted to discussion about the obtained 
results, and the last section provides their summary. 

2 Analysed system and symmetry of rules 

We analyse a one-dimensional cellular automaton, with two possible states 
of each cell, and 1. The rule of the change of a state of a cell is determined 
on the basis of the state of the cell itself and its two nearest neighbours. 
Namely, the rule is a chain (A, B, ..., H) of outcomes of the states of these 
three cells, in a decreasing order (111, 110, ...,000). Such a definition leads 
to 2^ = 256 different rules of the evolution process, so called elementary 
cellular automata. It was however shown that the number of unique rules is 
lower, as some of the rules are equivalent [8] . The rules which are equivalent 
to a given one are obtained as its conjugation, reflection and combined oper- 
ation of conjugation and reflection (see Tab{T|. As a result, the set of all 256 
elementary automata is reduced into 88 unique rules, or in other words, to 
88 groups of equivalent rules. The same result was recently obtained in |9] 
by an analysis of the symmetry of the rules. It may happen however, that 
after any change, the obtained rule is the same as at the beginning. Because 
of that, symmetry groups do not always contain four rules. To be precise, 
there are 44 groups which contain 4 rules each, 36 groups which contain 2 
rules each, and a remaining 8 rules which form one-element groups. Illus- 
trative examples are presented in Tabj2j 

Consider a system consisting of all possible sequences of zeros and ones 
of a specified length N. Each sequence can be converted to an exact single 
sequence determined by the automaton rule which is currently being con- 
sidered. Each sequence may, however, be obtained by the transformation 
of different numbers of sequences. The exact number of states (sequences) 
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Table 2: Examples of groups of equivalent rules: a) Group contains 4 dif- 
ferent rules No 7,21,31 and 87, b) Group contains 2 different rules No 132 
and 222, c) Group contains rule No 77 

leading to a given state depends on the automaton rule and the system 
size A^. The state space considered can be represented as a directed graph, 
where 2 nodes are identified with the possible states of the system. As for 
a given automaton, each state can be transformed to only one other state 
the out-degree of all nodes is equal to 1. However, the in-degree is different 
for different nodes. 

3 Classes 

Let us consider a given automaton rule. Analysis of the pattern of ties in the 
obtained graph allows us to indicate nodes which are similar, which means 
that the structure of connections of some nodes is the same. Because of that, 
nodes can be divided into classes, and each class contains nodes which are 
similar. As a result, the state space formed by all possible classes is smaller 
than the initial network of states. The rate of the reduction depends on the 
system size. 



The concept of the reduction of the system size is based on the observa- 
tion that, often in real systems, some kind of symmetry can be uncovered. If 



it is the case, it is possible to specify groups of states of the analysed system 
which manifest similar properties. When one thinks about the system like 
a graph, the mentioned similarities concern primarily nodes degrees. How- 
ever, this quantity by itself does not reflect the node properties. If two nodes 
have the same degree, but the degrees of their neighbours are different, the 
nodes cannot be classified as being similar. If weights of particular links in 
a graph are not equal, which is not the case here, one must also distinguish 
identical patterns of connections but with different weights. 

The class is defined as the set of nodes which have the same number of 
neighbours which belong to the same classes. The method of classes (groups 
of similar nodes) identification was described in detail in our previous pa- 
pers [ini in [H]. However, in the case of cellular automata the procedure is 
slightly different from that previously used. The underlying reason for the 
difference is that now the graph of transitions between states is directed. 
Therefore, the class of states is determined both by the list of states which 
for a particular automaton lead to a given state, and the one-element list, 
which preserve information about the state which is obtained from the one 
actually being considered. The parameter which differentiate two nodes is 
their degree, so, at the very beginning of our procedure, all nodes which 
have the same degree are labelled with the same symbol. With this done, 
the next step is to check if the lists of symbols of neighbours of all nodes 
which have assigned a given symbol are exactly the same. If not, appropri- 
ate elements have to be distinguished: for each class where this ambiguity 
is observed, as many separate classes must be introduced as it is necessary 
to ensure that for any two states, which belong to a given class, the full 
agreement for both lists - states which lead to a given state and the state 
which are obtained from it - is obtained. In general, the obtained network 
of classes is weighted, as more than one of the neighbours of a given node 
in the network of states may belong to the same class. In a further analysis 
of the obtained network of classes one must take into account the number 
of states which belong to particular classes. The procedure usually leads to 
a significant reduction of the size of the system. 

It must be emphasised, that for each system size A'^, the class structure 
obtained for the rules known to be equivalent [51 19] should be exactly the 
same. An interesting question is if it is also true for the rules which belong to 
different symmetry groups, but equivalent according to the function ^(N). 



4 Results 

4.1 Number of classes in a function of system size for par- 
ticular automata 

The number of classes for particular automata changes with the system 
size N. We can indicate a few different types of this dependence: 

• automata No and 255 

the number of classes is constant regardless of the system size and 
equals 2; one class contains states composed of zeros or ones, respec- 
tively, and the second class contains all remaining possible sequences. 
All of them are transformed in one case to the state composed of zeros, 
and in the other to the state composed of ones. In other words, in this 
case we obtain the star graph 

• automata No 15, 51, 85, 170, 204 and 240 

the second characteristic group of rules contains those for which just 
one class is obtained, irrespective of the system size; the rules cause 
that both in and out degree of each state is equal to 1; for the rules 
which belong to this group a new state of each cell is obtained as 
follows: for rule No 204 a new state is equal to the previous one, for 
rule No 51 it is negation of the previous state; for the four remaining 
rules a new state is determined on the basis of the previous state of the 
one of the nearest neighbours, and it is the state currently assigned 
to the appropriate neighbour or its negation, namely it is: for rule No 
170 the state of the right neighbour, for rule No 85 the negation of the 
right neighbour, for rule No 240 the state of the left neighbour, and 
for rule No 15 the negation of the left neighbour 

• automata No 105 and 150 

the number of classes in the case of these two rules form a sequence 
presented in TabjSJ during evolution in accordance with rule No 105, 
if the number of ones in the triple of cells (a given cell and its two 
nearest neighbours) is even - or 2 - irrespective of the previous state 
of a given cell, it is set to 1, if it is odd the value is set to 0. In the 
case of rule No 150 the situation is opposite - the new value is equal 
to if the number of ones is even, and to 1 if it is odd 

• automata No 154, 166, 180 and 210 

in this case, if N is an odd number, the number of classes equals 1, 
whereas if it is even the number of classes form a sequence presented 
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Table 3: Number of classes found for the automata No 105 and 150 in a 
function of the system size A^. 
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Table 4: Number of classes found for the automata No 154, 166, 180 and 210 
in a function of its even A^, in N in odd number of classes equals 1. 



in TabJ4| it may be noted that for every second even value (indivisible 
by 4) the number of classes is equal to 5; for all four rules classified 
here we observe the four configurations of a given cell state and its 
neighbours lead to the change of the cell state, and in particular if 
all three states are the same, the state of the considered cell does not 
change 

automata No 45, 75, 89 and 101 

here, similarly to the previous case, if A^ is odd the number of classes 
is equal to 1, but if N is even the number of classes increases expo- 
nentially with N for A^ ^ 10; the rules classified to this category are 
a mirror image of the rules described above. Because of that, obser- 
vations about the number of changed states remain the same, but in 
this case, if states of all cells in the triple are the same, the state of 
the cell in the middle is flipped 

automata No 60, 102, 153 and 195 

for those four automata if the system size is an odd number, the num- 
ber of classes equals 2, and if the system size is even the number of 
classes changes with A^ in the way presented in TabjSj for every second 
even value (indivisible by 4) the number of classes equals 3; all rules 
classified here are symmetric with respect to its center, and also in this 
case four configurations of states in the triple of cells lead to change 
and to maintain the current cell state 

automata No 90 and 165 

the pattern of change of the number of classes in the case of those two 

automata is particularly interesting, more because the rule produce in 
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Table 5: Number of classes found for automata No 60, 102, 153 and 195 in 
a function of its even N, in N in odd number of classes equals 2. 
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Figure 1: Number of classes in a function of the system size for automata 
No 90 and 165. 



the evolution the Sierpinski triangle fractal [T] ; the number of classes 
obtained in this case is presented in Fig{T| in this case, half of the 
triples also lead to the modification of the current cell value, and, 
for rule No 90, a new value is equal to 1 if the value of one of the 
nearest neighbours has the same value as the currently considered cell, 
otherwise it is equal to 0; for rule No 165 the situation is opposite 

remaining automata 

for most automata, none specified above, the number of classes in- 
creases exponentially with the system size, although slower than the 
number of possible sequences. The number of classes grows at the 
most as approximately 1.7^, while the number of states as 2^, where 
A^ is the system size. The largest number of classes is observed for 
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Figure 2: Histogram of the number of classes n for A^ = 19 in the case of 
the group of the automata for which exponential increase of the number of 
classes with the system size is observed. Note that 2^^ = 524288, then the 
reduction is by at least one order of magnitude. 



the automaton No 30, which belongs to the third class in the Wolfram 
classification system, and displays chaotic behaviour |Jj . Of course the 
same number of classes was identified for 3 remaining automata which 
belong to the same symmetry group (in a Wolfram's sense), it is the 
automata No 86, 135 and 149. The second largest number of classes 
is observed for the automata No 110, 124, 137 and 193, known to be 
Turing complete [U [13] . Those automata demonstrate a mixed char- 
acter of evolution, where regular and irregular behaviour is observed. 
The reduction of the size of the system for different automata is very 
different, and changes from one to even four orders of magnitude. 
The obtained values should be compared with 2 , i.e. the number of 
states. A histogram of the number of classes n for A^ = 19 is presented 
in Figj2} 

An exponential increase of the number of classes with the system size 
is very often observed in the whole range of analysed lengths, and only 
in some cases a deviation from this rule appears for automata shorter 



than approximately 10 cells. Among automata for which an exponen- 
tial increase of the number of classes is observed, a group for which 
the number of classes depends on the system size parity is noticeable. 

All but two rules for which atypical (different that exponential) character 
of the change of the number of classes with the sequence size are classified 
in j5] as surjective rules. The rules which do not fulfil this property are 
automata No and 255. What's more, most of them are also classified as 
being chaotic in accordance with the definition of chaos presented in pL4j . 
For all those automata value of z-parameter U^ is equal to 1, which is 
expected as they are chaotic. However, our classification indicate that for 
two symmetry groups in the Wolfram sense which are also classified as being 
chaotic, namely the groups [30, 86, 135, 149] and [106, 120, 169, 225], a typical 
- exponential - increase of the number of classes with the sequence size is 
observed. 

4.2 Networks of states and classes 

Our results show that the same number of classes can be found for automata 
which belong to different groups, out of 88 possibilities; in most cases this 
happens for a small number of groups for small N . Whether the structure 
of classes is the same or the observed compliance concerns only the number 
of classes, depends on the system size. However, in some cases the equiva- 
lence seems to be independent of A^. The detailed description is presented 



in Tab, 11 The first column contains 88 groups of rules, one after another. 
The following columns present the numbers of classes identified, as depen- 
dent on the system size N . The number of classes is supplemented with a 
letter, which makes it possible to distinguish cases for which the number of 
classes is the same but the lists of neighbouring states (in- and out-states) 
are different. The same class structure for different groups is observed or 
not, depending on N . 

For automata for which the same symbol was assigned, the complete 
consistency of the networks of states and the (reduced) networks of classes 
is observed. As a criterion of the equivalence of networks of states we take 



the identity of lists of in- and out-neighbours of all states. In Tab, 11 it is 
seen, for example, that the symbol lA for N = 5 appears 8 times, while 
for A = 18 only 4 times. The result means that in a previous case for 
8 symmetry groups, namely (15,85), (45,75,89,101), (51), (105), (150), 
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(154,166,180,210), (170,240) and (204), networks of states and classes for 
A^ = 5 is exactly the same. The same is true for A^ = 18 for groups (15, 85), 
(51), (170,240) and (204). The case of automata No 51, 204 and groups 
(170,240) and (15,85) should be considered separately. In this case, both 
networks - states and classes - are the same (1^) regardless of the system size 
N. This is because these automata are reversible: the identity and negation, 
the left- shift and the right-shift and their negation. There, states take part 
in long (e.g. left-shift) or short (identity) loops. In principle, there is no 
difference between states. Then, each state belongs to the same class, and 
there is only one class. Hence, each state is produced by a state from the 
same class and also gives a state from the same class. Hence, the reduced 
network is exactly the same. 

On the other hand, it is worth noting that in the case of a few groups, 
the pattern of the change of the number of classes with the system size is 
exactly the same. This is true in the case of automata No 23 and 232, 77 
and 178, 105 and 150, the pair (43, 113) and the pair (142,212), the group 
(12, 68, 207, 221) and the group (34, 48, 187, 243). A conjecture appears, and 
this equivalence is an indication of some symmetry of the rules, different 
than the conjugation and the reflection. Apart from operations which allow 
us to indicate equivalent rules, one can also consider negation and rever- 
sion of rules TabJH} Those operations lead to sequences, which belong to 
the same symmetry group, usually different than the groups to which the 
original sequence belongs. For 16 elementary automata, the sequence ob- 
tained as a negation and as a reversion is exactly the same. Those automata 
form 8 one-element and 4 two-elements symmetry groups TabJT} The group 
which contain the transformed sequence is joined with the group which con- 
tain the original one. We also observe that two four-elements groups which 
contain sequences obtained in accordance with the negation and reversion 
operations, are joined (Tabjsj) because of the same pattern of change of the 
number of classes. 

If we take into account the pattern of change of the number of classes, 
we see that the number of symmetry groups is equal to 80 (the number of 



different rows in Tab, 11). However, it is possible that this equivalence is 



broken for larger systems. 
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Table 6: Rule symmetry. The bar means negation, i.e. = 1 and 1 = 0. 
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Table 7: Equivalence of the selected symmetry groups. Single horizontal 
lines separate different symmetry groups, and double lines separate groups 
which are joined. 
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Table 8: Equivalence of the selected symmetry groups. The single horizontal 
line separates different symmetry groups. 
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Table 9: Number of groups as dependent on the system size A'^. 

4.3 Number of groups in a function of system size 

It is obvious that the structure of our classes is exactly the same for all 
rules which belong to each from among 88 unique symmetry groups [S] |9] . 
However, our class identification procedure allows us to take into account, 
not only the symmetry which characterizes the automata rules [H [9], but 
also the symmetry which manifests in the structure of transitions between 
particular possible sequences of a given length. Here, a group consists of 
automata for which the number of classes #{N) and the network of states 
as dependent on the system size N is the same. This is the basis for our 
classification. As a result, the number of groups for any system size is less 
than 88. The number of groups in a function of the system size is presented 
in Tabi] 

4.4 Cycles 

If the number of classes is equal to 1, both in and out degree of each node 
equals 1. It means that there is a sequence of configurations in which, start- 
ing from a given configuration, passing through other possible configurations 
of the system, we are able to return to the initial state. We examine the 
length of such cycles for the automata for which one class was obtained. As 
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Table 10: Lengths of cycles found for the automata No 45 on lengths A^ = 15 
(number of classes equals 1). Here, Lc is a cycle length, iic is the number of 
cycles of a given length and m shows how many times A^ is housed completely 
in Lc 

a result, histograms of the cycle lengths for different automaton lengths were 
constructed. In most cases the length of the cycle is an integer multiple of 
the system size. 



An example is presented in Tab, 10, where the result of the number of 
cycles of a given length for the rule No 45 is shown for A^ = 15. Symbols 
used in the table specify: Lc is a cycle length, Uc is the number of cycles of 
a given length and m shows how many times N is housed completely in Lc- 
The example presented is representative and it shows all possible situations 
which may be observed. First of all we can observe the situation where 
the length of the cycle is very short, in the example presented it is 1,2 or 
3 configurations. The length of the cycle equal to 1 means that the given 
configuration passes into itself. Then, we can observe a few cases where the 
length of the cycle is a multiple of the configuration length, in the case of 
the automaton No 45 for A^ = 15, it is the case in 6 cases. In most cases 
we also observe one cycle which appears as many times as the length of the 
configuration does (in Tab 10 marked as *). Finally, lengths of cycles which 



cannot be easily interpreted may also arise (in Tab, 10 marked as o). As we 
have seen above, a cycle often is equivalent to one class. 



5 Discussion 

We present a new method of cellular automata classification which allows us 
to differentiate rules in respect of the character of connections between all 
possible states of the system of a given size. We show that different types 
of character of the change of the number of classes with the system size is 
observed for the elementary CA. For most of them, an exponential increase 
is observed, but we also observe some number of rules which show different 
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kinds of behaviour. Namely, besides the obvious case of automata No 
and 255, a non-trivial dependence of the number of classes in a function 
of the system size is observed for automata No 15,45,60,90,105 and 154 
(and their equivalent rules). In most cases, the class identification allows 
a significant reduction of the system size. The proposed method indicates 
that the number of symmetry groups is equal to 80, as for some groups 
indicated by Wolfram, the same pattern of the change of the number of 
classes for different system size is observed. The proposed method also 
makes it possible to show that number of symmetry groups changes with 
the system size. 

The class reduction speeds up the calculation of the stationary probability 
distribution of particular configurations [12 . In particular, all „Garden of 
Eden" configurations have the probability zero, and all configurations which 
form a cycle have the same probabilities. In our method, it appears that all 
states which belong to a given class have the same probability. Hence, once 
we are able to calculate the probability of a given class we also know the 
probability of each state which belongs to this class. As a consequence, any 
analysis which concerns stationary properties of the analysed system may be 
performed on a reduced, smaller representation of the system. The method 
was described in pjj. Although the results presented here are obtained 
for the elementary cellular automata, our method can be applied to any 
family of automata, irrespective of its rules and size of its neighbourhood. 
It should be also emphasised that the method is general, and it can be 
applied to any system for which we are able to define a state space and an 
elementary process which transforms one state into another. As our method 
is numerical and not analytical, one should be cautious in predicting an 
asymptotic behaviour of the system. What we are able to conclude, is the 
behaviour of the system in the range of finite N. 
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